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Frequency  Response  ef  Multipass  SheH-and-Tube  Heat  Exchangers 


LEWIS  ISCOL  R.  I  ALTPETER 


NOMENCLATURE 

The  following  nomenclature  Is  used  in  the  paper: 


a, 

a^ 


=  ^  (dl) 


hg  A 

yigCf 


(dl) 


A  =  heat  transfer  area  per  pass  (sq  ft) 


b,  = 


h^  A 


(dl) 


''K 


(dl) 


,  V. 

=  heat  capacity  (BTU/lb  F* ) 


C,  C'  =  arbitrary  constants  (dl) 
d  =  thermal  diffusivity  of  pipe  wall  (sq  ft/hr) 
D  =  intermediate  parameter  (dl) 
e  =  2.7182818.,.  (dl) 


E  =  Intermediate  parameter  (**F) 
f  =  intermediate  parameter  (dl) 

F  =  cross  sectional  area  of  metal  in  pipe  wall  (sq  ft) 
g  =  intermediate  parameter  (dl) 
h  =  heat  transfer  coefficient  (BTU/hr  sq  ft  P*) 
k  «=  thennal  conductivity  (BTU/hr  ft  P*) 

K  *  Intermediate  parameter  (dl) 

L  =  tube  length  per  pass  (ft) 
m  =  number  of  shell  passes  (dl) 
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M  «  Intermediate  parameter  (dl) 

n  ■  number  of  tube  passes  in  each  direction  per  shell  pass  (dl) 
N  =s  intermediate  parameter  (dl) 
p  «  intermediate  parameter  (dl) 

P  »  heat  transfer  perimeter  of  pipe  (ft) 
q  »  constant  of  integration  (®p) 


.21-  (dl) 


R 

s 

S 

t 


intermediate  parameter  (dl) 
Laplace  transform  variable  (dl) 
shellside  fluid  temperature  (•?) 
Laplace  transform  of  S  (•?) 
time  (hr) 


T  =  tubeside  fluid  temperature  (*p) 

T  =  Laplace  transform  of  T  (•?) 

U  =»  overall  heat  transfer  coefficient 
V  »  fluid  velocity  {ft/hv) 

V,  V  =  definite  integrals  (*p) 
w  =  fluid  flow  rate  (Ib/hr) 

W  =  pipe  wall  temperature  (®p) 

=  Laplace  transform  of  W  (®P) 

X  =  axial  distance  along  exchanger  (ft 
y  =  axial  distance  along  exchanger  (dl 
Z  =  intermediate  parameter  (dl) 


U  A 


w,  c^ 


U  A 

c. 


(dl) 


(dl) 


(BTU/hr  sq  ft  F*) 
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•^r  =  pipe  wall  thickness  (ft) 

S  =  Intermediate  parameter  (dl' 
i  =  root  of  auxiliary  cubic  (dl) 

©  =  time  (dl) 

A  =  radial  distance  through  pipe  wall  (ft) 

?  =  Intermediate  parameter  (dl) 

=  density  (ib/cu  ft) 

=  radial  distance  through  pipe  wall  (dl) 
f  fo*  =  dummy  variables 
V'  -  Intermediate  parameter  (dl) 

Subscripts; 

1  =  ordinal  iiumber  of  shell  pass 

j  =  ordinal  number  of  tube  pass  within  a  particular  shell  pass 
Q  =  summation  Index 
R  =  reference 
S  =  shellslde 
T  =  tubeside 
W  =  pipe  wall 

As  usual  the  use  of  an  Independent  variable  as  a  subscript  of 
a  dependent  variable  Indicates  partial  differentiation, 

as  = 
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iT 

Tx 


I 

INTRODUCTION 

A  multipass  shell  and  tube  exchanger  may  be  described 
by  a  set  of  partial  differential  equations  similar,  in  many 
respects,  to  those  describing  a  counterflow^''’*^^  or  one  side 
lumpecf^^  exchanger.  This  set  of  equations  may  be  thought  of 
as  constituting  a  mathematical  model  of  the  exchanger. 

Models  of  varying  degrees  of  complexity  may  be  constructed 
for  the  same  exchanger.  In  this  paper  models  are  first 
constructed  which  neglect  the  heat  capacity  of  tube  walls. 
Refinements  are  then  made  which  allow  the  walls  to  be 
introduced  as  either  lumped  or  distributed  thermal  capacity. 

The  models  presented  herein  are  quite  formidable  of  aspect, 
involving  the  definition  of  many  sets  of  Intermediate  parameters. 
It  should  be  borne  in  mind  that  the  question  which  must  be 
answered  regarding  the  feasibility  of  use  of  a  given  model 
is  not,  "How  complicated  is  the  model?"  but,  "How  much  does 
it  cost  to  extract  the  desired  information  from  the  model?". 

The  desired  information  is  here  the  frequency  response 
characteristics  of  the  exchanger.  This  Information  may  be 
particularly  easily  extracted  if  an  explicit  expression  is 
obtainable  for  the  transfer  function.  Such  an  expression 
is  obtainable  for  the  models  considered.  In  fact,  it  is  the 
object  of  this  paper  to  show  how  such  explicit  transfer 
functions  may  be  obtained. 

The  evaluation  at  a  particular  frequency  of  a  transfer 
function  of  the  complexity  here  considered  takes  on  the  order 
of  15  seconds  on  an  Intermediate  speed  digital  computer  and  so 

^  Nximbers  in  parentheses  designate  References  at  end  of  paper. 
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Ptg.l  Fluid  flow  Inside  pipe. 
Direction  of  flow  determines  sign 
of  Tx  term  In  equation. 


K - ' 

ORIENTATION  2 

Fig. 2  The  two  possible  orientations  of 
a  1-2  exchanger. 


costs  about  iO.25.  This  cost  would  be  cut  at  least  50  per  cent 
by  the  use  of  a  high  speed  computer.  Clearly,  then,  the  use 
of  these  models  even  for  routine  work  is  definitely  feasible. 

FUNDAMENTAL  EQUATIONS 

The  partial  differential  equation  describing  the  temper¬ 
ature  changes  of  a  fluid  flowing  inside  a  pipe  (see  figure  la) 
la''"  : 

'1'  ^  7  (S  -  T) 

If  the  direction  of  flow  of  the  fluid  is  reversed  (figure  lb) 
this  becomes: 

(2)  -  T,  (3.T) 

The  assumptions  embodied  in  these  equations  are  essentially 
those  employed  to  develop  the  usual  expression  for  the 
log  mean  temperature  dlfference'^^  and  are  as  follows: 

T  '  » 


1.  U,  w, 


and  V  are  constant 


2.  Plug  flow  prevails, 

3*  No  temperature  gradients  exist  in  the  fluid  other 
than  in  the  axial  direction. 

4.  No  partial  phase  changes  take  place  in  the  system. 
Equations  (1)  and  (2)  are  made  dimensionless  by  the  introduction 
of  new  space  and  time  variables, 

(5)  y=^ 

(4) 6  =1^ 

L 

Using  (3)  and  (4),  (l)  becomes: 

(5)  +  —  Te  =  -  T) 

Each  pass  of  a  multipass  exchanger  is  described  by  an  equation 
similar  to  (5).  These  equations  together  with  appropriate 
boundary  conditions  constitute  one  model  of  the  exchanger. 

THE  1  -  2n  EXCHANGER 

4 

Two  orientations  of  a  1-2  exchanger  are  possible  (see 
figure  2) .  For  each  orientation  there  are  four  possible 
temperature  forcing  transfer  functions.  Forcing  may  be 
applied  on  either  the  tube  or  shell  side,  and  the  response 
may  be  taken  on  either  side. 

Transfer  functions  will  be  derived  for  the  case  of 
orientation  1,  forcing  on  the  tube  side,  with  response  on 
either  side.  Disturbances  frequently  take  the  path,  forcing 
on  the  tube  side,  response  on  the  tube  side.  Corrective 
action,  however,  frequently  Involves  flow  forcing  which  is 
not  here  considered.  The  other  six  possible  temperature 


forcing  transfer  functions  may  be  derived  by  procedures  very 
similar  to  the  one5  to  be  demonstrated. 

The  equations  for  the  l-2n  exchanger  (see  figure  5) 
are,  in  dimensionless  form: 


(6) 

(Tjjj  + 

(Ti)s 

=  ‘^.(s 

-  T-) 

J  =  1,3,5,. . .2n-l 

(7) 

(Ti). 

=  ‘“.(S  - 

Ti) 

J  =  2,4,6, . . .2n 

(8) 

Sj  +  r 

s® 

=  ^ 

(Ti  -  S) 

Note  that  tubeside  velocity  is  chosen  as  reference  velocity. 
The  additional  assumption  is  made  that  the  area  of  each  pass 

is  the  same.  Boundary  conditions  are: 


(9) 

Tj  (t,0)  =  T-., 

(t,0) 

J  =  3,5,7 

(10) 

Tj  (t,l)  =  T-., 

(t,l) 

j  =  2,4,6 

(11) 

S(t,0)  =  0 

The 

solution  is  carried  out  in  detail 

in  Appendix  A. 

THE  GENERAL  MULTIPASS  SHELL  AND  TUBE  EXCHANGER  (m  -  2mn) 

Under  the  additional  mild  restriction  that  there  be  no 
heat  transfer  between  different  shell  passes  (also  a  usual 

(if) 

steady  state  assumption  )  it  is  possible  to  derive  explicit 
transfer  functions  for  the  general  multipass  shell  and  tube 
exchanger  having  m  shell  passes  and  2n  tube  passes  per. shell 
pass.  The  notation  is  illustrated  by  figure  4.  When  two 
subscripts  are  used  the  first  always  refers  to  the  number  of 
the  shell  pass.  The  equations  take  the  form: 
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y-0 


Fl£.3  Heat-exchanger  model  with  lumped 
thermal  capacity  In  wall. 


Pig. 4  A  3-12  exchanger  Illustrating 
the  notation  used  for  mltlpass  ex¬ 
changers.  Many  other  orientations  are 
possible,  but  all  are  similarly  bandied. 


(12)  +  (Ti-  +  (Tij)e  =  ^»(Si  -  Ti]  ) 

(13)  +  (Si  )^  +  r  (Si)^ 

i 

Boundary  conditions  mirror  the  continuity  of  fluid 
before.  Details  of  the  method  are  in  Appendix  B. 


flow,  as 


LUMPED  THERMAL  CAPACITY  IN  WALLS 

Thus  far  models  have  been  considered  that  have  no 
thermal  capacity  in  the  pipe  walls.  Now  it  will  be  shown 
how  lumped  thermal  capacity  may  be  introduced  into  the  pipe 
wall  between  any  two  fluids.  Furthermore  it  will  be  shown 
that  the  equation  describing  the  more  complicated  system 
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may  be  reduced  to  the  form  previously  considered. 

The  system  of  figure  5  been  described  by  the  equations 

(14)  (3  -  T) 

(15)  +  r  -  S) 

It  may  also  be  described  by  the  set: 

(16)  Ty  +  Tq  =  a.  (W  -  T) 

(17)  +  r  Se  =  a,,  (W  -  S) 

(18)  W©  =  b.  (T  -  W)  +  bj,  (S  -  W) 

Transforming  and  rearranging  (18)  one  obtains: 

(19)  W  = - T  +  - ? 

s  +  b,  +  b;^  s  +  b,  +  b^ 

Using  this  expression  for  W  in  the ’transforms  of  (I6)  and  (17) » 

a,  b,  n  _  a,  b,. 

s  +  a, - - - T  =  - - - - -  S’ 

s  +  b,  +  bj^  s  +  b,  +  b^ 

(21)  +  ^rs  +  a^  - 

Equations  (20)  and  (21)  are  of  the  same  form  as  the  transforms 
of  (l4)  and  (15),  but  with  different  constants.  The  solution 
from  this  point  on  is  the  same  as  before. 

The  treatment  for  the  case  of  a  wall  in  contact  v.’ith 
a  single  fluid  (as  a  shell  wall)  is  similar.  It  should 
be  noted  that  in  a  single  exchanger  the  the^’nial  capacity  of 


some  walls  may  be  neglected,  the  thermal  capnclty  oC  other 
walls  lumped,  and  the  capacity  of  still  other  walls  treated 
In  distributed  fashion  as  will  be  described  below. 


DISTRIBUTED  THERMAL  CAPACITY  IN  WALLS 

The  system  of  figure  6  has  been  described  by  the  equations; 


(22)  ^  (S  -  T) 


(23)  Sy  +  -1  Sf  =  JLE -  (T  -  S) 

V?  wjr 

It  may  also  be  described  by  the  set: 


(24)  Tx  +  -1  Tf  =  _ \w(+,x,0)-tV 

w.,  c,  L  -1 

(25)  S,  =-^^|^W(t,x,^)  -  s\ 


Pl^.6  Heat-exchanger  model  vltb  dis¬ 
tributed  thermal  capacity  In  wall. 


(26)  =  d 

The  boundary  conditions  for  (26)  are; 

(27)  k  A  (W^  A  ^W(t,x,0)  -  T^ 

(28)  k  A  (V)^  ;^_  =  hj.  A^  S  -  W(t,x,^)] 

The  independent  variables  are  made  dimensionless  by  the 
transfoinnations  (3),  (4),  and, 

(00)  o-  =  _ 

The  reference  velocity  is  v.^  .  A^’ter  employing  the  Laplace 
trensforrr.at  j  on  and  reari^anging,  the  set  takes  the  form: 

(30)  Tv^  +  (s  +  OC.  )  T  =  W(s,y,0) 

(31)  +  (rs  +  )  3  -  ‘»‘,.V(s,y,l) 

(3?)  Iw  -  =  0 
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(33) 


V,jg(s,y,0) 


^3T  + 


(34)  ('C  U 

Equation  (52)  may  be  solved  using  boundary  conditions  (33) 
and  (34)  to  obtain  W(s,y,o-  )  as  a  linear  combination  of 
S,  W(s,y,0),  and  l^(s,y,l).  If  this  expression  for  W(s,y,. 
is  evaluated  at  o-  =  o  and  =1,  one  obtains: 

^35)  ’  W(c,y,0)  +  ^i.~(s,y,l)  =  '^3?  +  5,^3 

(36)  ’^«W(s,y,0)  +  'i^W(s,y-,l)  = 

where , 

(37)  r,  =  ^  (e*  -  e*  )  +  V,  (e* 

(38)  .  0 

(39)  =  V,  (e 

(40)  S,  =  9 

(41)  =  2», 

(4a)  -«-®)  + 

(43)  1^3  =  aV, 

(44)  >4,=  y^(e'°  +  e-® 


+  e 


) 


iS  -y?  . 

+  ^  ) 


+  e' 


) 


Prom  the  linear  set  (35)  and  (36)  one  obtains  a  solution 
of  the  form: 

(45)  VJ^(s,y,0)  =  D,  T  +  Dx  S 

(46)  W(s,y,l)  =  Dg  T  +  S 


Substitution  from  (45)  and  (46)  into  (30)  and  (31)  reduces 
the  system  of  equations  to  the  form  previously  ’considered. 

Walls  vihich  are  in  contact  with  a  single  fluid  may  be 
similarly  handled.  Any  number  of  walls  in  an  exchanger  may 
be  handled  in  distributed  fashion  as  shown  here. 

The  reader  will  have  noted  that  the  Fourier  heat  conduction 
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equation  is,  used  in  rectangular  coordinates  rather  than 
cylindr:’ cal ,  as  is  strictly  required.  -The  error  introduced 
is  not  large  and  is  outweighed  by  the  convenience  of  being 
able  to  obtain  an  explicit  transfer  function  in  terms  of 
the  elementary  functions. 

It  should  be  noted  that  models  in  the  literature  of 
counterflow^*'*"^  ,  parallel  flow^''*'^  ,  and  one  side  lumpecF'^ 

exchangers  may  be  refined  and  thus  undoubtedly  brought  into 
better  agreement  with  experimental  work  by  the  addition  of 
distributed  thermal  capacity  in  the  pipe  walls,  as  provided 
for  here . 

APPENDIX  A 

Equations  (6),  (?)»  and  (8)  may  be  transformed  and 
put  into  the  form: 

(*^7)  (^i  )y  +  Tj  =  g,  S  ^  I,  j  =  1,3,5, ...2n-l 

(48)  -(T-  )y  +  f,  =  g,  S  -  J  =  2,4,6,...2n 

(49)  ^  Ti 

where,  ' 

(50)  f,  =  s  + 

(51)  ft  =  2no‘-t  +  ns 

(52)  g.  = 

(53)  g*.  = 

The  boundary  conditions  become: 


(54) 

fj  (s,0)  = 

(s,0) 

J  -  3,5,7,** 

(55) 

Tj  (s,l)  = 

Tj.»  (s,l) 

j  =  2,4,6,.. 

(56) 

S(s,0)  =  0 

Equations  (4?)  and  (48)  possess  the  well-known  solutions 
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(57)  Tj 


«3  —  1  f  •  •  • 


2n-l 


=  e 


f. 

^  <3y  +  Qj  e 
-f'l 


-fj 


,  .  A,  c  ' 

(58)  Tj  =  J  g,  5  dy  +  q;  e'*’’  J  =  2,4,...2n 


The  Qj  may  be  evaluated  In  order.  For  instance  directly 
from  (57): 

(59)  q.  =T,  (s,0) 

Then  since, 

(60)  (s,l)  =  T.  (s,l) 

may  be  obtained,  etc.  It  may  be  verified  that: 


1-P 

(61)  q.  =  V  -t; -  +  V 

'  ^  l_l  -  E_J 


mi-p 


1-P 


3 


jA 


+  p  ^  J  =  1,5» • • .2n-l 


(62) 


+  pK  J  =  2,4, 


.  .2n 


where 

(65) 

(64) 

(65) 


T  (s,0)  has  been  set  to  unity  and. 


For  the  extension  of  this  method  of  derivation  to  the  general 
multipass  case  it  will  be  Important  that  the  qj  are  linear 
combinations  of  V,  v' ,  and  1  with  coefficients  of  the  form 
^  ,  where  the  C^.  and  Cq  are  constants. 

The  7:  may  be  eliminated  from  (49)  using  (57)  and  (58) 
and  the  resulting  equation  differentiated  twice  to  remove 
the  Integrals.  One  obtains: 
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(6m)  -  C  +  (?r,  nf,  -  f*-  )S  =  0 

Three  Initial  conditions  are  needed.  ^{s,0)  is  obtainable  from 
(56).  S,^  (s,0)  and  (s,0)  are  obtainable  by  setting  y  =  0 

before  the  first  and  second  differentiations,  respectively^^  . 

(s,0)  and  (s,0)  turn  out  to  be  linear  combinations 

of  the  and  thus  linear  combinations  of  V  ,  V*  ,  and  1  with 

coefficients  of  the  same  form  as  before.  The  algebraic  details 
will  be  presented  at  the  end  of  the  anpendix  so  as  not  to  confuse 
the  development. 

The  solution  of  (66)  is  of  the  form: 

(6?)  S  =  E,  e^*^  +  e*^-‘i  +  eM 

where  are  the  roots  of  the  auxiliary 

cubic.  These  roots  may  most  simply  be  found  using  Cardan's 
formula''^^  .  Ej  ,  ,  and  E3  are  related  to  the  initial 

conditions  on  by  the  linear  set; 


> 

The  E  are  thus  also  linear  combinations  of  V  ,  V  ,  and  unity. 
One  may  define: 


(69)  E 

\  = 
'L 

K,  V  + 

X 

M  , 

X 

+ 

N. 

3 

3 

3 

s 

.Eauation 

(67)  may  now  be  rewritten; 

('"0)  S’  = 

K. 

h-k. 

e’‘^ 

+ 

K3 

+  v'j^ 

1  ^ 

e**"i 

+ 

M3 

-  1 

+ 

N, 
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^  ,  and  ^'^are  defined  from  (70)  by: 


(71)  s  =  )^,  +  5^v+ 

The  are  also  of  the  form  Cq  e^'l  .  All  that  must 

now  be  done  before  'S’(s,l)  (the  shellside  response)  or  (s,0) 
(the  tubeside  response)  can  be  evaluated  numerically  is  to 
evaluate  V  and  V*  .  This  may  be  accomplished  using  (71). 

The  symbol  (  ^  »  defined  by: 

(72)  (  /»  ,  ^o')  =  r  /• 

Equation  (71)  is  multiplied  by  g  e'*^^  and  integrated  with 

t 

respect  to  y  between  0  and  1.  Using  the  notation  defined 
by  (72)  this  yields: 


(73)  V  =  (  ^.  ,  g^  )  +  V(J^  ,  g,  )  +  v'  (  g,  e  '^'^  ) 

fyj 

Similarly,  if  (71)  is  multiplied  by  g^  e  '  : 

(74)  v'  =  (y&.,  gj  e'^^  )  +  v(yi  ,  g,  e^'i  )  +V*()^',  g,  e'^^  ) 

The  integrations  indicated  in  (75)  and  (74)  may  be  easily 
carried  out  due  to  the  simple  form  of  the  .  The  second 
order  set  of  linear  algebraic  equations  defined  by  (75) 
and  (74)  may  then  be  solved  numerically  for  V  and  v' 

With  V  and  v'  in  hand  ^(s,l)  is  available  numerically  from  (71) 
and  (s,0)  from  (58)  and  (62). 

All  that  remains  is  to  present  expressions  for  the  ^ 
as  linear  combinations  of  V  ,  ,  and  unity;  and  to  show  the 

results  of  the  integrations  indicated  in  (73)  and  (7^)* 

This  will  be  done  using  several  sets  of  Intermediate  parameters. 


Z 
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Z*  ,  ,  and  Z^  are  similarly  defined  with  (-f,  )  substituted 

for  (f.  ). 

/  \  O  1  “  P 


(76)  R.  = 

(77)  Ri.  = 

(78)  R,  = 


1  -  p 

Pd-P'^  ) 

1  -  P 

n  1  -  p'^ 

1  -  p"  (1-p)^ 


n  P(1  -  P'"  ). 

(79)  R.  =  - 

’  1  -  p  (1  -  p)*- 

Finally,  the  K,  M,  and  N  defined  by  (69)  are  given  by: 


(80)  =  Z,  R3  +  z;  Ry 

i  z  2 

3  3  3 


(81)  M,  =  Z.  PR3  +  z;  pR^ 


(82)  N.  =  Z,  R.  +  Z’.  Rv 

Z  Z  2 

3  3  3 

The  are  now  given  by  (70)  and  (71).  The  integrals  follow. 


■  -1 


Ki-. 

e-^3 

)  and  (  ® 

) 

are 

similarly  evaluated  using 

K. 

,  K3  and  M,  , 

Mx 

> 

M3  ,  respectively,  for 

N.  . 

Nr  . 

• 

,  e^'^  )  are  similarly 

(y. 

),  ('f> ,  ), 

and 

(Z'* 

evaluated  with  (f,  )  substituted  for  (-f^  )  . 


APPENDIX  B 

Equations  (12)  and  (15)  may  be  transformed  and  put  into 
the  form: 


(84)  +(T..  )  +  r  -  g, 


(85)  +(?;  )>.  +  3,  = 


T- 


Equations  (84)  possess  known  solutions 


(f) 


odd  j  i 
even  [ 

f  odd  r] 

1  even  1 


(86) 


=  + 


S'i  dy  +  q.- 


The  q^j  may  again  be  evaluated  in  order.  They  will  be 
found  to  be  linear  combinations  of  2m  definite  integrals 


(V,  , 

V.’  , 

.  .  •  •  •  .  Vi  ) 

f-  -sa 

I  S-  e  dy 

and  unity. 

(87) 

V. 

\ 

J 

(88) 

/ 

V/ 

4 

=  J 

* '  fu 

g;  e  •  J  dy 

from 

► 

(86)  may  be  substituted 

into  (Be)  and  the  resulting 

equations  differentiated  tv/ice  to  eliminate  the  integrals.  This 
yields,  for  each  of  the  ,  a  third  order  differential  equation. 

The  solution  may  be  written: 


(89) 


+ 


E 


tZ. 


+ 


The  E;  are  not,  however,  linear  combinations  of  only 
(V.  ....  V*  ,1).  Because  of  the  boundary  condition 

(90)  S’,-  (s,0)  =  (s,0)  (odd  i  ,  1  1) 

(91)  (s,l)  =  s^_,  (s,l)  (even  1), 

the  Et  are  linear  combinations  also  of  (s,0)  or  (s»l) 
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depending  on  whether  i  Is  odd  or  even  respectively.  Rearranging 
(89)  as  was  done  to  (6?)  to  obtain  (71): 

(92)  5.  =  yT.,  V,‘  +  V,  +  V,'  +  . . 

+  (s,?) 

The  equation  for  has  no  term  In  ^  (s,0).  In  the 

equation  for  the  term  In  {s,l)  may  be  eliminated  by 

substitution  from  the  equation  for  S,  .  Proceeding  stepwise 
the  following  set  is  obtained  from  (92). 

(93)  =  Ko+'%l^y;i  V;  +  V.'J 

Each  of  the  m  equations  of  the  set  (95)  may  be  multiplied 

fey  S,  and  integrated;  and  then  by  g,  e**'!  and  integrated 

as  was  done  to  (71),  to  obtain  the  set  of  2m  equations: 


(94)  V.  =  ,  g,  e'^'i) 


e  ) 


+  v; 


j-i 

f" 

(95)  V.'  =  (J^,-^  ,  g,  1  V;  (Y':j.g,  )  +  V.'  (J^-/^  , 

This  set  may  be  solved  for  the  2m  constants  (V^  ,  ...  } . 

The  desired  responses  are  then  easily  evaluated,  as  before. 


e 
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